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Abstract 

We derive the exact equations of motion (in Newtonian, F — ma, 
form) for test masses in Schwarzschild and Gullstrand-Painleve coordi- 
nates. These equations of motion are simpler than the usual geodesic 
equations obtained from Christoffel tensors in that the afhne parameter is 
eliminated. The various terms can be compared against tests of gravity. 
In force form, gravity can be interpreted as resulting from a flux of su- 
perluminal particles (gravitons). We show that the first order relativistic 
correction to Newton's gravity results from a two graviton interaction. 

1 Gullstrand-Painleve Coordinates 

In general relativity, the Schwarzschild solution for a spherically symmetric (non 
rotating) black hole has been known since 1915. The usual choice of coordinates 
is the one Karl Schwarzschild used in its discovery, Schwarzschild coordinates, 
which are characterized as keeping the metric diagonal, but have a coordinate 
singularity at r = 2M: 

(dr) 2 = (1 - 2M/r) dt 2 - dr 2 /(l - 2M/r) - r 2 {d9 2 + sm 2 (6) d</> 2 ), (1) 

where we have chosen coordinates with G = c = 1. Note that if we multiply r 
and t by M, the metric will end up with an overall multiple of M 2 which we 
can cancel. For convenience, we will do this both with the Schwarzschild and 
GP coordinates. The reader can reinsert M by making the reverse substitution. 

Gullstrand-Painleve (GP) coordinates were discovered by Allvar Gullstrand3 
[J and Paul Painleve in 1921/1922: 

dr 2 = (1 - 2M/r)dt 2 - 2y/2M/r dt dr - dr 2 - r 2 (d8 2 + sm 2 (0) d(/) 2 ). (2) 



"This essay received an "honorable mention" in the 2009 Essay Competition of the Gravity 
Research Foundation. 

^^Gullstrand had a primary role in denying Einstein a Nobel prize for relativity. 
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While the metric is not diagonal, the curvature is concentrated into the dt 2 and 
dr dt terms. The purely spatial terms of the metric, — dr 2 — r 2 (d0 2 + sin 2 (9) d(j) 2 ), 
are identical to the spatial part of the natural metric for a flat (Minkowski) space. 
This makes GP coordinates a natural choice for a model of the gravitational force 
by an exchanged particle (which we will call the graviton). 

A major step towards unifying general relativity with quantum mechanics 
was made by Anthony Lasenby, Chris Doran, and Stephen Gull in 1993-1998 
when they invented "Gauge Theory Gravity" (GTG). [3] In their 

theory, they rewrote the tensor theory of general relativity into the language 
of Dirac's gamma matrices [5] [H] (also called, with various subtle shades of 
meaning, the "Spacetime Algebra" , "Geometric Algebra" , or "Clifford Alge- 
bra"). This enables calculations of interactions of fermions with black holes. 

The GTG differed from general relativity only in that it was built on a flat 
(Minkowski) background space and consequently could not support topologi- 
cally interesting solutions to Einstein's field equations. Those who might doubt 
the existence of science fiction topics such as wormholes would find the GTG a 
substantial improvement over standard general relativity. 

The application of GP coordinates to the GTG is not obvious from their 
papers; a short explanation may be useful. The metric for a black hole in GTG, 
requires the definition of four ancillary functions of radius, /i, f 2 , <?i, and g 2 
and (see equation (52) of Ref. [3j ) is: 

dr 2 = (1 - 2M/r)dt 2 - 2(f l92 - f 29l )dr dt-(f 2 - f 2 )dr 2 - r 2 (d6 2 + sin 2 (6) dcj> 2 ). 

(3) 

The four functions of radius define the directional derivatives with respect to 
r and t. One has a choice of gauge for the radial direction that allows one to 
choose g 2 arbitrarily, and from this compute /i, f 2 , and g\. One requires that 
/i and gi go to one at infinity, while f 2 and g 2 approach zero. To obtain GP 
coordinates, one would use the freedom in g 2 to require that /f — /f = 1. This 
gives a metric whose spatial portion is flat and therefore is GP. See Refs. [12] 
and [13] for generalizations of GP to the Kerr metric. 



2 The Force of Relativistic Gravity 

When computing orbits in general relativity, the easiest and traditional method 
is to use Christoffel tensors T^ A . This gives four differential equations in the 
affine parameter q: 

dq 2 + vX dq dq ~ W 

for fj, = 0, 1, 2, and 3. For massive particles, the affine parameter q can be taken 
to be the proper time but this fails in the massless case. Newtonian gravitation 
avoids the use of q and so is able to get by with just three differential equations 
in t: 

dV _ d$ n 
m ^ + Gm d X J=°> (5) 
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that is, ma — F = 0. 

To find the Newtonian equations of motion from the Schwarzschild metric, 
the most direct way is to first note that geodesic paths extremize s = J ds. So 
one can write ds in terms of dx/dt, dy/dt, dz/dt and use the Euler-Lagrange 
equations to vary J ds. Following the Newtonian tradition, we will abbreviate 
dx/dt by x and similarly for x, etc. The resul1@ for Schwarzschild coordinates 
is a few terms in r and r — 2: 

x r 4 (r — 2) = 'ixr 2 r 2 — x(r — 2) 2 + 4yr(yx — xy) + 4zr(zx — xz) + , 
2r(r — 2)y(yx — xy) + 2r(r — 2)z[zx — xz), 

and similarly for y and z. GP coordinates are simpler in that they need only 
powers of r: 

xr 5 = —xr 2 + 2xr — 2r 2 y{xy — yx) — 2r 2 z(xz — zx) + 3xr 2 r 2 + , 
\/2r(3r 2 x — r 3 x\r\ 2 + l.5r 3 xf 2 + 2y(xy — yx) + 2z(xz — zx)) 

The easiest way to verify the above equations is to choose a set of random 
positions and velocities, and compare the acceleration with that computed from 
the geodesic equations. One finds that the above are exact, and that they work 
for massless as well as massive test particles. 

We can modify general relativity by making changes to the various terms 
in the above. This gives us information on which terms have been fixed by 
experimentaal tests. For example, the terms contribute to the small deflection 
of starlight (DoS) and perihelion advance of Mercury (PoM) in proportion as 
follows: 



GP term 


DoS 


PoM 


— x/r [i 


+ 1/2 





2x/r A 





-1/3 


-2y(xy - yx)/r 3 


+1 


+4/3 


3xr 2 /r 3 


-1/2 





3V2x/r 2 - 5 








-V2±|rf /r 1 - 5 








l.bV^xP/r 1 - 5 








\/8y(xy - yx)/r 4 - 5 








Total: 


+1 


+1 



(8) 



Only the first four terms contribute to the solar system tests of general relativity. 

A characteristic of a spherically symmetric gravitational force is that it de- 
pends only on the distance to the body r, and two velocities, the raidal velocity 
r and the velocity perpendicular h. For the Schwarzschild coordinates, we find: 

r = (-l-2h 2 )/r 2 + 3r 2 /(r(r-2)) + 2/r 3 , 

h = 2rh/r 2 +4rh/{r 2 (r -2)). ( ) 

2 The author used the symbolic computation software MAXIMA. A Java applet that 
demonstrates various properties of relativistic gravity that uses these formulas is at 
www.brannenworks.com/Gravity 
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The — 1/r 2 is the Newtonian part of the force. Other than the +2/r 3 , the 
remaining terms are all of 2nd order in j3 = y 1 x 2 + y 2 . 

For GP coordinates, grouping terms by their order in r we find: 

r = -V2r/i 2 /r 1 ' 5 + (-l-2/i 2 + 3r 2 )/r 2 + 3V2r/r 2 - 5 + 2/r 3 , 
h = -V^hyr 1 - 5 + 2rh/r 2 + V2h/r 2 - 5 + ^/T/2hr 2 /r 4 - 5 . 

For r = h = 0, Schwarzschild and Painleve coordinates give the same "static" 
acceleration: 

r = -l/r 2 + 2/r 3 , 

h = 0. 1 J 

This acceleration differs from the Newtonian acceleration by the addition of the 
+2/r 3 term. In particular, the static acceleration is zero at the event horizon, is 
negative inside the event horizon, and is everywhere identical for both GP and 
Schwarzschild coordinates. 



3 Gravitons in the Gravitostatic Limit 

The classical electric force between charged bodies follows an exact 1/r 2 law. 
This can be attributed to the electric force being carried by an elementary 
particle (the photon) that is massless and whose intensity naturally decreases 
according to the area over which the photons are spread. In making these sorts 
of arguments, we must recognize that we cannot assume that the velocity of 
the gravitons involved is the same as the speed of light. In particular, in GP 
coordinates, a particle falling through the event horizon exceeds the speed of 
light (that is, \r\ > 1). 

Gravitons capable of producing such a force must also exceed the speed of 
light. Such a theory would have gravity waves also travel faster than light. The 
experimental measurement of the speed of gravity is a subject much discussed 
in the physics literature. The mainstream view is that the speed has not yet 
been measured. [T51 Consequently, our discussion, while speculative, is not yet 
eliminated by experiment. 

Let the number of gravitons passing through a sphere of radius r be N(r). 
From Eq. (|TT|) . we have that 

N(r) cx |47rr 2 (-l/r 2 + 2/r 3 )| = 4tt(1 - 2/r). (12) 

As the radius increases, N(r) increases until it approaches a constant corre- 
sponding to Newtonian gravitation. We will use n for this constant so that 
N(r) = (1 — 2/r)n. The number of gravitons per unit volume at a distance r 
from the gravitating body is 

p(r) = N{r)/{Anr 2 ) = (1/r 2 - 2/r 3 )n/(47r). (13) 

For large r, this approaches n/iAnr 2 ). 
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The number of gravitons appearing at a radius of r is given by the derivative 
of N(r): 

N'(r) = 2n/r 2 . (14) 



Dividing this number by the surface area over which these gravitons appear, 
4irr 2 , gives the rate of change of gravitons per unit volume: 



The number of gravitons appearing at the radius r, Eq. (|15[) . is proportional to 
1/r 4 . This is approximately proportional to the square of the graviton density 
at that radius, Eq. (Ti"3"|) . Consequently, we obtain that the graviton density 
increases at a rate proportional to the square of the density of gravitons. In a 
perturbation theory of gravitons, this term would arise by allowing two gravitons 
to interact to create three (or more) all moving in the same direction. That is, 
in addition to their emission from massive bodies, gravitons are emitted by 
stimulated emission somewhat similar to the effect in lasers. 
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